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ABSTRACT 

We use our new stellar population models, which include effects from variable element 
abundance ratios, to model the Ca4227 absorption line indices of early-type galaxies 
(Trager et al. 1998), and to derive calcium element abundances. We find that calcium, 
although being an cv-element, is depressed with respect to the other a-elements by up 
to a factor 2. This confirms quantitatively earlier speculations that early- type galaxies 
are calcium under abundant. We find a clear correlation between a/Ca ratio and central 
velocity dispersion, which implies that more massive galaxies are more calcium under- 
abundant. Interestingly this correlation extends down to the dwarf spheroidal galaxies 
of the Local Group for which a/Ca ratios have been measured from high-resolution 
spectroscopy of individual stars (Shetrone et al.). The increase of the calcium under- 
abundance with galaxy mass balances the higher total metallicities of more massive 
galaxies, so that calcium abundance in early-type galaxies is fairly constant and in 
particular does not increase with increasing galaxy mass. This result may be the key 
to understand why the Call triplet absorption of early-type galaxies at 8600 A is con- 
stant to within 5 per cent over a large range of velocity dispersions (Saglia et al. 2002; 
Cenarro et al. 2003). The origin of the calcium underabundance in early- type galaxies 
remains yet to be understood. We argue that formation timcscalcs are disfavoured 
to produce calcium underabundance, and that the option of metallicity dependent 
supernova yields may be the most promising track to follow. 
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evolution - galaxies: formation - galaxies: stellar content 



1 INTRODUCTION 

In the late seventies it has first been noticed that early-type 
galaxies show relatively weak Call triplet (CaT) absorption 
at 8600 A (Cohen 1979; Faber & French 1980). The metal- 
licities derived from these lines lie well below solar, which 
seems hard to reconcile with the super-solar metallicities 
obtained from the classical Lick indices Mg 2 and (Fe) (e.g., 
Trager et al. 2000a; Kuntschner et al. 2001; Thomas et al. 
2002; Maraston et al. 2003). Moreover, the very tight corre- 
lation between Mg 2 and velocity dispersion a (e.g., Bender, 
Burstein & Faber 1993) has not been found for the CaT line 
(Cohen 1979; Faber & French 1980). Subsequent measure- 
ments of CaT absorption in early-type galaxies providing 
increasingly better data quality confirm these early results 
(Terlevich, Diaz & Terlevich 1990; Peletier et al. 1999). Very 
recently, Saglia et al. (2002) measured the CaT absorption 
in 94 early- type galaxies at very high signal-to-noise ratios, 
using an improved index definition (Cenarro et al. 2001). 
They again find a relatively weak CaT absorption and even 



a mild anti-correlation between CaT and a with a remark- 
ably small scatter (see also Cenarro et al. 2003). 

Saglia et al. (2002) show that neither uncertainties in 
the CaT modelling nor an (already previously suggested) 
steepening of the initial mass function provide compelling 
explanations. A very attractive way out remains to be the 
conclusion that early-type galaxies are indeed underabun- 
dant in the element calcium. The main problem is that this 
interpretation remains speculative as long as the sensitivity 
of the CaT index to calcium abundance is not known, and in 
particular calcium underabundant stellar population models 
of CaT are not available. 

Interestingly, also the blue Lick index Ca4227 at A w 
4227 A is very weak in early-type galaxies, which again in- 
dicates a possible calcium underabundance (Worthey 1992; 
Vazdekis et al. 1997; Trager et al. 1998). A similar conclu- 
sion has been drawn on the basis of Ca4455 index measure- 
ments (Worthey 1998; Henry & Worthey 1999). We note, 
however, that Ca4455, despite its name, appears to be only 
sensitive to the abundances of iron-peak elements (Tripicco 
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& Bell 1995), so that the weak Ca4455 absorption reflects 
the well-known iron underabundance in early-type galaxies 
(e.g., Worthey, Faber & Gonzalez 1992). Ca4227, instead, is 
sensitive to calcium abundance, and — different from CaT — 
for Ca4227 we now have a set of stellar population models at 
hand (Thomas, Maraston & Bender 2003, hereafter TMB) 
that explicitly take abundance effects of individual elements 
into account, and in particular are computed for different 
calcium abundances. Here we use these models to analyse 
the data of Trager et al. (1998), and for the first time to 
quantitatively derive calcium element abundances of early- 
type galaxies. 

In Section 2 we briefly introduce the stellar population 
model. The results are presented in Section 3 and discussed 
in Section 4. 



2 MODELS 

In TMB we present stellar population models for the first 
time with different chemical mixtures and element abun- 
dance ratios. We computed all optical Lick indices from CNi 
to TiO-2 in the wavelength range 4000 < A < 6500 A of Sim- 
ple Stellar Population (SSP) models varying the abundance 
ratios of a-elements (i.e. O, Mg, Ca, Na, Si, Ti) to iron peak 
elements (Fe, Cr) and to the individual elements carbon, 
nitrogen, and calcium. The impact from the element abun- 
dance changes on the Lick absorption-line indices are taken 
from Tripicco & Bell (1995). The models comprise ages be- 
tween 1 and 15 Gyr and metallicities between 1/200 and 5 
solar. We refer to TMB for more details. 

Particular care was taken to calibrate the SSP mod- 
els on globular cluster data, which — most importantly — 
include objects with relatively high, namely solar, metal- 
licities (Puzia et al. 2002). We match very well their Ca4227 
indices with models in which calcium is enhanced like the 
other a-elements with respect to the iron-peak elements, 
in agreement with results from high-resolution spectroscopy 
of individual stars (see TMB). Note that calcium itself is 
an a-element. Having the issue of weak calcium absorp- 
tion in early-type galaxies in mind, we developed addi- 
tional models in which calcium is detached from the group 
of a-elements, allowing for variations of the a/Ca element 
ratio. In TMB, models of the index Ca4227 for the ele- 
ment ratios [a/Ca] = —0.1, 0.0, 0.2, 0.5 are provided, with 
[a/Ca] ee log {X a /X c .) - log (X a /X Ca ) Q . 

We will use these models to derive a/Ca ratios of early- 
type galaxies. Note that the above notation relates the ele- 
ment ratio to the solar value. We will speak of calcium under- 
abundance, if a/Ca is larger than solar, hence [a/Ca] > 0. 



3 RESULTS 

We analyse the early-type galaxy data of Trager et al. 
(1998). This turned out to be the only data set in the liter- 
ature with Ca4227 indices for 'normal' early-type galaxies. 
We do not use the sample of Longhetti et al. (2000), because 
it consists of peculiar, disturbed and shell galaxies. Still, it 
should be kept in mind that the Trager et al. (1998) sample 
contains a rather inhomogeneous morphological mix, with 
many SO and later-type galaxies included. From this sample 



we selected those objects, for which the error in the Ca4227 
measurement is smaller than 0.25 A, which yields 39 galaxies 
with a median error of 0.2 A. 



3.1 [MgFe]' vs. Ca4227 

In Fig. 1 we show Ca4227 as a function of the index [MgFe]', 
with [MgFe]' ee ^Mg b ■ (0.72 • Fe5270 + 0.28 • Fe5335) as 
defined by TMB. [MgFe]' is independent of a/Fe ratio vari- 
ations, and therefore a good tracer of total metallicity. SSP 
models for the a/Ca ratios [a/Ca] = —0.1, 0.0, 0.2, 0.5 and 
the metallicities [Z/E] = -0.33, 0.0, 0.35, 0.67 at the fixed 
age t = 12 Gyr are plotted in the left-hand panel (see labels 
in the figure). Models of constant metallicity and varying 
a/Ca ratio (dashed lines) are almost perfectly vertical, in- 
dicating that [MgFe]' is not sensitive to calcium abundance. 
Solid lines in the left-hand panel of Fig. 1 are models with 
constant a/Ca ratios. Note that the models in the [MgFe]'- 
Ca4227 plane are highly degenerate in age, as both indices 
have similar age dependencies. This is shown with the right- 
hand panel of Fig. 1, in which models with ages 1, 4, and 
15 Gyr and metallicities [Z/R] = -0.33, 0.0, 0.35, 0.67 are 
plotted (see labels in the figure), the a/Ca ratio being fixed 
at the solar value. The left-hand panel of Fig. 1 therefore is 
well suited for isolating a/Ca ratio effects. 

The galaxy data are shown as filled circles, the open 
circle is their median. The filled triangle is the integrated 
light of the Galactic Bulge (Puzia et al. 2002) . The data lie 
below the model with solar a/Ca ratio, no combination of 
age and metallicity can be found to match the data. a/Ca 
ratios from solar to roughly 3 solar, < [a/Ca] < 0.5, are 
obtained. The maximal calcium underabundance relative to 
the other a-elements thus is about a factor 2-3. We are 
confident that this detection of calcium underabundance is 
significant in spite of the relatively large observational er- 
rors for two reasons. 1) The Ca4227 range covered by the 
data plotted here is consistent with the data of Longhetti 
et al. (2000) , which have, although being disturbed and shell 
galaxies, smaller observational errors (< 0.09 A). 2) The me- 
dian of the sample (open circle) is at Ca4227 = 1.1±0.03 A, 
significantly below the model with solar a/Ca. From the lo- 
cation of the median we infer that calcium is depleted in 
early-type galaxies by typically a factor 1.4. 

The large triangle in Fig. 1 is the average of 15 Galac- 
tic Bulge fields measured by Puzia et al. (2002). Obviously, 
the Galactic Bulge is also calcium underabundant. Its a/Ca 
ratio is consistent with the median of the early-type galaxy 
sample, at somewhat lower total metallicity. This result is 
in good agreement with the detailed abundance analysis of 
individual field Bulge stars from high-resolution spectra by 
Mc William & Rich (1994). These authors measure super- 
solar a/Ca ratios ([a/Ca] ~ 0.2), hence calcium underabun- 
dance, for the metal-rich stars ([Fe/H] > —0.5) of their sam- 
ple. 

The galaxies M31 and M32 are labelled in Fig. 1. M31 
has stronger Ca4227 and [MgFe]' indices than M32, so that 
both objects have very similar a/Ca ratios ([a/Ca] ~ 0). 
Curiously, M32 shows stronger CaT absorption than M31 
(Kormendy & Bender 1999), which seems to indicate higher 
calcium abundance in M32. Note however, that the CaT 
index is severely contaminated by Paschen absorption lines 
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Figure 1. Early-type galaxies in the [MgFe]'-Ca4227 plane. Galaxy data (filled circles) are from Tragcr et al. (1998). The open circle 
is their median. The filled triangle is the integrated Galactic Bulge light from Puzia et al. (2002). Lines are stellar population models 
from Thomas et al. (2003). Left-hand panel: Models for constant a/Ca ratios [«/Ca] = —0.1, 0.0, 0.2, 0.5 (solid lines) and constant 
mctallicities [Z/H] = —0.33, 0.0, 0.35, 0.67 (dashed lines) at a fixed age t = 12 Gyr. Right-hand panel: Models for constant ages of 1, 
4, and 15 Gyr (solid lines), and the constant metallicities [Z/H] = —0.33, 0.0, 0.35, 0.67 (dashed lines) at a fixed solar a/Ca clement 
ratio ([a/Ca] = 0). 



(Cenarro et al. 2001). As M32 is younger than M31 (Trager 
et al. 2000b), the strong CaT absorption index of M32 is 
most likely caused by an increase of Paschen line absorption. 



3.2 Correlations with galaxy velocity dispersion 

In Fig. 2, the index Ca4227 (top panel) and the a/Ca ratios 
derived from Fig. 1 (bottom panel) are plotted as functions 
of central velocity dispersion (circles) . More precisely, we de- 
termine with the TMB models ages, metallicities, and a/Ca 
ratios from the indices H/3, [MgFe]', and Ca4227. It should 
be emphasized, that uncertainties in the stellar population 
ages have only a minor effect on the determination of the 
a/Ca ratio as demonstrated by Fig. 1. Velocity dispersions 
are taken from Whitmore, McElroy & Tonry (1985), Faber 
et al. (1989), Bender et al. (1992), and Gonzalez (1993). The 
triangles in Fig. 2 are the Local Group dwarf spheroidal 
galaxies Carina, Draco, Fornax, Leo I, Sextans, and Ursa 
Minor. Their mean a/Ca ratios are the average of the el- 
ement abundances of individual stars (~ 5 per object) de- 
termined from high-resolution spectroscopy (Shetrone et al. 
2001; Shetrone et al. 2003). These authors derive a-element 
abundances by averaging the abundances of the elements 
Mg, Ca, and Ti. Central velocity dispersions are taken from 
the compilation of Mateo (1998). 

Ca4227 correlates with velocity dispersion within the 
measurement errors. The increase of Ca4227 with increas- 
ing <7o, however, is balanced by the simultaneous increase 
of [MgFe]', so that also a/Ca increases with velocity disper- 



sion. A linear least-square fit (solid line) to the giant galaxy 
data yields 

[a/Ca] = -0.03 + 0.0008 * a . (1) 

In other words, more massive galaxies are more metal-rich 
and more calcium underabundant. The admittedly large 
scatter of the relation between a/Ca and oo is fully consis- 
tent with the errors in [a/Ca] (~ 0.1 dex), which are caused 
mainly by the large observational errors of Ca4227. Ca4227 
is a rather weak absorption line, and requires high signal-to- 
noise to be measured accurately. Given these large errors, it 
seems even surprising that such a clear correlation could be 
found. Most interestingly, it is very well consistent with the 
a/Ca ratios of the Local Group dwarf spheroidals. We note 
that the large spread in Fig. 2 may also come from the in- 
homogeneous morphological mix in the Trager et al. (1998) 
sample. Ca4227 index measurements for a similarly large 
sample of elliptical galaxies at significantly higher signal-to- 
noise than the Trager et al. (1998) data would be very useful 
to check the present result. 

The increase of [a/Ca] by ~ 0.2 dex between ao = 100 
and 300 km/s is comparable to the typical metallicity in- 
crease (Trager et al. 2000b; Kuntschner et al. 2001, D. 
Thomas et al., in preparation) in that mass range. Hence 
calcium abundance stays approximately constant or may 
even decrease (the present data quality does not allow us 
to discriminate these possibilities). The Ca4227 index in- 
creases, mainly because of its sensitivity not only to cal- 
cium abundance but also to total metallicity (Tripicco & 
Bell 1995). The Paschen line corrected Call triplet index 
CaT* at 8600 A (Cenarro et al. 2001), instead, is found to 
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Figure 2. Ca4227 index (top panel) and a/Ca ratio (bottom 
panel) of early-type galaxies as functions of central velocity dis- 
persion. Ca4227 data are from Trager et al. (1998, circles). a/Ca 
ratios are derived from Fig. 1. Triangles are Local Group dwarf 
spheroidal galaxies for which a/Ca ratios are measured in indi- 
vidual stars (Shctrone et al. 2001; Shctronc ct al. 2003). Central 
velocity dispersions are from Whitmore et al. (1985), Faber et al. 
(1989), Bender, Burstcin & Faber (1992), Gonzalez (1993), and 
Mateo (1998). The solid line is a linear least-square fit to the giant 
galaxy data (circles). 



be constant to within 5 per cent over the same range of ve- 
locity dispersions (Saglia et al. 2002; Cenarro et al. 2003). 
This result can be understood now, provided that the index 
CaT* is more sensitive to calcium abundance than Ca4227. 



populations (Matteucci 1994; Thomas, Greggio & Bender 
1999). Note that high a/Fe ratio are actually an iron under- 
abundance (Buzzoni et al. 1994; Trager et al. 2000a, TMB). 
They are a good measure for star formation timescales, 
because a substantial fraction of iron comes from the de- 
layed enrichment by Type la supernovae (Greggio & Renzini 
1983; Matteucci & Greggio 1986; Thomas, Greggio & Bender 
1998). Note that also the old, metal-poor stars in the halo 
of the Milky Way, which formed at an early stage of the 
Galaxy's evolution, are underabundant in iron (McWilliam 
1997, and references therein) . It may seem therefore straight- 
forward to link also the underabundance of calcium with 
formation timescales. 

The situation for the a/Ca ratio, however, is entirely 
different. The old, metal-poor halo stars are not underabun- 
dant in calcium (McWilliam 1997), hence the solar a/Ca 
ratio is already at place in the very beginning of chemical 
enrichment. In other words, calcium, which is actually an a- 
element, is enriched in lockstep with the other a-elements. 
Moreover, Type la supernovae produce mainly iron and do 
not contribute significantly to the enrichment of calcium 
(Nomoto, Thielemann & Yokoi 1984). The calcium under- 
abundance found in early-type galaxies can therefore not 
be explained through the link of formation timescales with 
delayed Type la supernova enrichment. 

Within Type II supernova nucleosynthesis, calcium 
comes mainly from supernovae of the less massive progen- 
itors with m < 20 M© (Woosley & Weaver 1995; Thiele- 
mann, Nomoto & Hashimoto 1996). If the chemical enrich- 
ment is dominated by stars more massive than ~ 20 M Q , 
calcium underabundance can be achieved (Molla & Garcfa- 
Vargas 2000). Because of the short lifetimes of these stars, 
however, extremely short formation timescales (~ 10 7 yr) 
would be required, implying unreasonably high star forma- 
tion rates up to 10 5 Mq per year. 

Alternatively, also a flattening of the initial mass func- 
tion (IMF) enhances the contribution from high-mass stars 
to the chemical enrichment. Integrating the supernova yields 
of Woosley & Weaver (1995) with variable IMF slopes, we 
find that an abundance ratio [a/Ca] ~ 0.2 requires ex- 
tremely shallow IMF slopes of the order x ~ 0.2 (Salpeter 
being x = 1.35). A flattening of the initial mass function to 
produce calcium underabundance can therefore be ruled out, 
as such shallow IMFs violate other observational constraints, 
in particular predict stellar mass-to-light ratios (Maraston 
1998) well above the dynamical determinations (Gerhard 
et al. 2001). 



4 DISCUSSION AND CONCLUSIONS 

In the previous section we show that early-type galaxies are 
calcium underabundant. The reason for this chemical pecu- 
liarity is not yet understood. In this section we will briefly 
discuss possible mechanisms that can lead to the formation 
of calcium underabundant stellar populations. 

4.1 Formation timescales and IMF 

More massive early- type galaxies have higher a/Fe ratios 
and older average ages (e.g., Trager et al. 2000b; Thomas 
et al. 2002; Terlevich & Forbes 2002; Caldwell et al. 2003), 
which indicates faster formation histories of their stellar 



4.2 Metallicity-dependent yields 

A more appealing solution may be connected to a metallicity 
dependence of supernova yields, in the sense that in metal- 
rich supernovae the production (or at least the ejection) of 
calcium is depressed (Worthey 1998). 

From the theoretical side of Type II supernova mod- 
elling, there are no clear indications for the existence of 
metallicity-dependent yields. The calculations of Woosley & 
Weaver (1995, model B) predict only a very mild increase of 
the a/Ca ratio with increasing metallicity in the supernova 
ejecta for the metallicity range 1/10,000 to solar. Extrapo- 
lating these yields to double-solar metallicity leads to [a/Ca] 
ratios, which are ~ 0.2 dex below the values derived here for 
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massive early-type galaxies. If metallicity effects are at work, 
they must set in at super-solar metallicities, a range which 
has not yet been explored in Type II supernova nucleosyn- 
thesis calculations. 

This option gets support from the fact that only the 
stars with metallicities above solar in the McWilliam & Rich 
(1994) data are calcium underabundant. It would also be 
very well consistent with the a/Ca-a correlation of early- 
type galaxies found in this paper. The most significant in- 
crease of the a/Ca ratio with galaxy mass occurs among 
giant galaxies in a regime of super-solar metallicity. The de- 
crease of the cv/Ca ratio when going to the Local Group 
dwarf spheroidals by at most 0.1 dex is relatively small, in- 
stead, given the huge drop in metallicity to about 1.5 dex 
below solar. To conclude, the possibility of metallicity de- 
pendent yields at super-solar metallicities remains a promis- 
ing track to follow in order to find the origin of the calcium 
underabundance in early-type galaxies. 
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